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ABSTRACT 

We introduce Coronal-Line Forest Active Galactic Nuclei (CLiF AGN), AGN 
which have a rich spectrum of forbidden high-ionization lines (FHILs, e.g. [FeVII], 
[FeX] and [NeV]), as well as relatively strong narrow (~300 km s“^) Ha emission 
when compared to the other Balmer transition lines. We find that the kinematics of 
the CLiF emitting region are similar to those of the forbidden low-ionization emission- 
line (FLIL) region. We compare emission line strengths of both FHILs and FLILs to 
CLOUDY photoionization results and find that the CLiF emitting region has higher 
densities (10^ ® < hh < lO^'® cm“^) when compared to the FLIL emitting region (10^ ° 
<nH < 10^-^ cm“^). We use the photoionization results to calculate the CLiF regions 
radial distances (0.04 < Kclif < 32.5 pc) and find that they are comparable to the 
dust grain sublimation distances (0.10 < Rsc/s < 4.3 pc). As a result we suggest that 
the inner torus wall is the most likely location of the CLiF region, and the unusual 
strength of the FHILs is due to a specific viewing angle giving a maximal view of the 
far wall of the torus without the continuum being revealed. 

Key words: galaxies:active — galaxies:Seyfert — quasarsiemission lines — 

quasarsigeneral 


1 INTRODUCTION 


In the standard active galactic nuclei (AGN) unification 
model ( Antonucci fc Miller||1985 I a large-scale height torus 
surrounds the inner AGN, obscuring our view of the central 
continuum source and the broad emission line region (BLR) 
from most directions. The inner wall of this torus lies at 
the sublimation radius of the most refractive dust - no dust 
can exist stably at smaller radii. This inner wall is heavily 
irradiated. It has been proposed that ablation of the sur¬ 


face produces a thermal wind of highly ionized gas (Krolik 
I&: Kriss|1995[ ). This gas may produce an unusual spectrum, 
which may include forbidden high ionization lines (FHILs), 
also known as ‘coronal lines’ l Appenzeller &: Wagner|]1991 1. 

From most orientations, either pole-on or edge on, the 
torus inner wall will be barely visible. However, there should 
be a small range of angles of which the torus only just hides 
the inner nucleus. For these type 2 (dust obscured) AGNs 
the rear inner wall will be maximally visible. Here we de¬ 
scribe a population of AGN that are candidates to be this 
special population - the Coronal Line Forest AGN (CLiF 
AGN). 


* E-mail:mrose@cfa.harvard.edu (MR) 


In this paper, we introduce the CLiF AGN as a distinct 
class of AGN, and present data for a sample of the first 
seven identified CLiF AGN. The cosmological parameters 
used throughout this paper assume a A cold dark matter 
(ACMD) cosmology with parameters Dm = 0.3, Ha = 0.7 
and Ho = 70 km s“^ Mpc“^. 


2 PROPERTIES OF CLIF AGN 


Most AGN show some weak (EW < 3 A) spectral lines from 
forbidden transitions of highly ionized ions (FHIL, ionization 
potential. Ip > 54.4 eV; the Hel l edge) in their spect ra, e.g. 
[NeV], [FeVII], [FeX], [FeXIV] ( [Penston et ZjligSlf] The 
presence of FHILs in the spectra of AGN unambiguously in¬ 
dicate the presence of an AGN, due to their high ionization 
potentials, reaching up to soft X-ray energies ~0.1-0.2 keV. 
FHILs are not generally associated with LINER nuclei or 
starbursts (Filippenko & Halpern 19841. FHILs are physi¬ 


cally distinct because their transitions have higher critical 


^ In this work we define Ip as the ionization potential required 
to achieve the current ionic stage by removing an electron from 


a lower-ionization stage I De Robertis fc Osterbrock|198^. 
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densities (nc ~ 10^'® cm Osterbrock & Ferland 200^ 


than lower ionization species (nc < 10 ® 


Osterbrock 


fc Ferland|[200^ . In most AGN the FHILs are blueshifted 


(Av ~-250 km s“ ) with respect to the rest frame of the 
AGN and have substantial velocity widths (FWHM ~1000 
km s“^), That are intermediate between those of the stan¬ 
dard low ionization forbidden lines (FLIL) of the narrow 
emission line region (NLR) (FWHM ~500 km s“^, e.g. 
[OIII]A5007), and the BLR (FWHM ~5000 km s"^). These 
kinematic properties support an origin for the FHIL emit¬ 


ting gas in a wind outflowing from the BLR (Appenzeller & 
Wagner|1991 1 or the torus I Krolik fc Kriss|1995[ ). 

A few rare AGN have instead a dozen or more FHILs 
with relatively high fluxes (F(FHIL)/F(H/3) > 0.25): we 
call these the CLiF AGN. So far, GLiF AGN have only 
been identified serendipitously. Indeed, there are just 5 GLiF 
AGN in the literature (Table [T|: Q1131-I-16 (Rose et al. 
Wlj ), HI Zw 77 ( |Osterbrock|[l98l|, M rk 1388 ( jOsterbrock 


(Fosbury & Sansom 


TMsI, ESQ 138 G1 (lAlloin et al.|199^ and Tololo 0109-383 


19831. 


The optical spectra of these objects are remarkable in 
several ways. As well as the the typical FLILs, i.e. [OHI], 
[Nil], [01] and [SH], and the Balmer emission lines including 
Ha, H/3, their spectra show many FHILs, including [FeV], 
[FeVI], [FeVH], [FeX], [FeXI], [NeV] and [ArV] all with 
F(FHIL)/F(H/3) > 0.2. This is more than found in typical 


AGN which just show the [FeVH]AA5721,60871ines (Appen- 
[zeller fc Wagner||199ip , and even these are weak (EW < 3 

AT 

We hnd ten unusual properties of CLiF AGN: 

1 . Detection of the [FeVH] lines 3759A, 5159A, 5720A 
and 6087A with the flux ratio F([FeVH]A6087)/F(H/3) > 
0.25. Where F(H/3) is the flux from the narrow H/3 emission 
line. 

2. F([FeX]A6374)/F(H/3) > 0.2. 

3. F([NeV]A3426)/F(H,d) > 1. 

4. F([FeV])/F(H/3) & F([FeVI])/F(H,S) > 0.2. 


5. Hq/H/ 3 > 2.9 [the Gase B value seen in the NLR; (Os- 
terbrock fc Ferland|2006| ] and H 7 /H/? ~ 0.47 (Case B), im¬ 


plying that the apparent level of dust extinction calculated 
using Ha/H/3 is greater than that calculated using H 7 /H/ 3 . 

6 . Higher level HI Balmer transitions (with wavelengths 
blueward of H5 and up to Hie; [Rose et al.|2011[ ) are present 
in emission. 

7. The FHILs are not blueshifted with respect to the 
FLILs (Av < 100 km s-^). 

8 . The velocity widths of the FHILs are narrow and the 
same as the FLILs (FWHM ~ 300 km s“^) using single 
Gaussian hts. 

9. The FHILs maintain their high EWs for longer than 3 
years. This is long compared with the observed emission line 
fading times for Type Hn supernovae ( Smith et al.|2009' I and 
tidal disruption events for stars falling on to a supermassive 
black hole (SMBH; Gezari et al.|200^ Komossa et al.|2009 


Wang et al. 20121, two astrophysical phenomena that also 


exhibit rich spectra of narrow emission lines including both 
FLILs and FHILs. 

10 . High EW (>3A) emission lines from un¬ 
known/unexpected species are present. 

This cluster of properties appear to be important for 
GLiF AGN. However, at this stage it is premature to for¬ 


mally define the GLiF AGN. We need a larger sample of 
GLiF AGN in order to produce a well-crafted definition. 


3 A SEARCH FOR CLIF AGN IN SDSS 


Given that CLiF AGN have very different spectral prop¬ 
erties compared to typical AGN, one might expect that CLiF 
AGN would be easy to select. We used each of the unusual 
properties listed above that could be selected from the SDSS 
DRIO databas43 

The properties of the emission lines in the SDSS spectra 


are determined using the BOSS pipeline (Bolton et al. 20121. 


Each line is modeled as a single Gaussian. The amplitudes, 
centroids, and widths of each emission line are optimized 
non-linearly to obtain a minimum- 7 ;^ fit to each emission 
line. All fitted emission lines are assumed to have the same 
redshift. The widths of the emission lines are assumed to 
have the same intrinsic velocity depending on the ‘width 
group’ they belong to. The relevant groups for the emis¬ 
sion lines used in this work are ‘Balmer’ (consisting of the 
Balmer emission/absorption features) and ‘emission’ (con¬ 
sisting of non-Balmer emission lines). The velocity widths 
are a strength weighted average over the width group. The 
continuum level is estimated either from the best-fit velocitw 
dispersion model if the spectrum is classified as a galaxjlH 
or from a linear £t either side of the emission line if the spec¬ 
trum is classified as a quasar. The parameters such as the 
fluxes, FWHM, redshifts, EW and continua levels, as well as 
their associated uncertainties, are all reported by the BOSS 
pipeline ( Bolton et al.|[20l2 K We use EW instead of a flux 
ratio with H/3 because our definition for a strong FHIL as the 
ratio of the FHIL flux strength to the flux from the narrow 
H/3 emission line is not usable with the SDSS database. As 
the SDSS database does not differentiate between BLR and 
NLR Balmer emission, we cannot perform searches based on 
the flux ratios. 

Six of the ten GLiF properties could not be searched 
using the BOSS database. Properties 1, 4 and 10 could not 
be selected because these Fe emission lines are not defined in 
the SDSS spectral fitting model (see Table 5 in [Bolt on et al.| 
(20121 for the full list of emission lines fitted by the SDSS). 


The kinematic properties 7 and 8 are not available in the 
BOSS database either. Also, there is no timing information 
to determine whether property 9 is observed. 

We attempted to select CLiF AGN candidates from the 
BOSS catalog, making use of the on-line query tool, for prop¬ 
erties 2, 3, 5 and 6 : 


• For property 2, Gelbord, Mullaney & Ward (20091 


noted that the [FeX]A6375 broad blue wing is often blended 
with the [OI]A6364 emission line which is included in the 
BOSS catalog. As [OIJ6300/6364 = 3 from atomic physics, 
if [FeX]A6375 is present and blended with [OI]A6364, it will 
reduce the [OI] 6300/6364 ratio. Hence we required that ob¬ 
jects had to have [OI] (6300/6364) < 3. A caveat is that 


^ http://skyserver.sdss3.org/drl0/en/tools/search/sql.aspx 
^ The best-fit velocity dispersion is determined by fitting locally 
for the position of the minimum-y^ versus trial velocity disper¬ 
sion in the neighborhood of the lowest gridded value. See 
https://www.sdss3.org/drl0/algorithms/redshifts.php for a full 
description. 
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Table 1. Basic properties of the CLiF AGN sample. Redshifts (z[q///]) are determined by single Gaussian fits to the prominent 
[OIIIJA5007 emission line. The column r’ gives the SDSS r’ model magnitudes (https://www.sdss3.org/drl0/algorithms/magnitudes.php). 
We choose model magnitudes because these objects are type 2 AGN and therefore the photometry will be dominated by the host galaxy. 
In the cases of ESO 138 G1 and Tololo 0109-383, there are no available r’ magnitudes, instead we use Gousins R band magnitudes taken 
from [Lauberts & Valentijn| ‘Spectrum’ gives the origin of the spectrum. The column ‘Ref’ refers to the first publication that 

discusses the optical spectrum of the CLiF AGN candidate. 


Name 

RA (J2000) 

Dec (J2000) 

^lOlll] 

r’ 

Spectrum 

Ref. 

Q1131-I-16 

11:31:11.05 

16:27:39.50 

0.1732 

17.31 

WHT/ISIS 

Rose et al. (2011) 

III Zw 77 

16:23:45.87 

41:04:56.69 

0.0341 

15.25 

SDSS 

Osterbrock (1981) 

Mrk 1388 

14:50:37.85 

22:44:03.61 

0.0216 

14.55 

SDSS 

Osterbrock (1985) 

ESO 138 G1 

16:51:20.13 

-59:14:05.20 

0.0091 

13.02“ 

ESO 

Alloin et al. (1992) 

Tololo 0109-383 

01:11:27.63 

-38:05:00.48 

0.0118 

12.38“ 

IPCS/RGO 

Fosbury & Sansom (1983) 

J1241+44 

12:41:34.25 

44:26:39.25 

0.0422 

15.98 

SDSS 

This work 

J1641-I-43 

16:41:26.91 

43:21:21.59 

0.2214 

18.60 

SDSS 

This work 


“ The R-band magnitude is from [Lauberts &: Valentijn] l |198^ . 


[OIJA6300 can also be blended with [SIIIJA6312. This blend 
conld result in a ratio [01] (6300/6364) > 3. However, the 
[SIII]A6312 emission line is included in the BOSS spectral 
fits ( [Bolton et al.|20l2 |. This selection criterion will then be 
conservative because any contribution from [SIIIJA6312 will 
result in CLiF AGN not being selected. 

• For property 3, we selected objects with strong 
[NeV]A3426 emission (E’V\f] > 10 A), although this only 
works for z > 0.109. 

• For property 5 we required that Ha/H/3 > 2.9 AND 
Hy/H/l < 0.5, so as not to exclude reddened objects. 

• For property 6 we note that He, which is in the BOSS 
spectral fits, has at A=3970 a wavelength close to that of 


[NeIII]A3968 even closer than [OIJA6364 and [FeX]A6375. 
Hence we created the criterion FW([NeIH]A3968) > 5 
as He will boost the [NeHI]A3968 flux. 


Each of the initial searches based on only one of these 
properties resulted in large samples of several thousand ob¬ 
jects. These were of unrelated phenomena, such as typical 
quasars whose broad emission lines altered both the Balmer 
decrements and [01] ratios, resulting in false detections for 
properties 2 and 5. Property 3 also selected quasars which 
lacked strong Fe FHILs. Normal type 2 AGN were also se¬ 
lected with property 3, suggesting that this emission line is 
important in type 2 AGN ( [Mignoli et al.|2013p . Finally prop¬ 
erty 6 selected both typical quasars and quiescent galaxies. 
These observations were confirmed with spot checks of their 
spectra. 

Searches based on pairs of the criteria also failed to con¬ 
strain the search. Therefore additional criteria were needed 
in order to constrain the search. 

We then searched for objects using all four properties 
(2 AND 3 AND 5 AND 6). This search produced just two 


* While EW has been used as a standard indicator of emission 
line strength, it depends on the intensity of both the emission 
line and underlying continuum. Therefore there is a caveat when 
defining the properties of CLiF AGN using the EW because the 
large observed EWs could be due to a weak stellar or AGN con¬ 
tinuum emission. 

® Tests of lower (~3A) EW limits returned larger samples of 
several thousand unrelated objects, and higher (~10A) EW limits 
produced no CLiF AGN at all. 


CLiF AGN. One of these was the previously detected Mrk 
1388, and one new CLiF AGN: J1641-I-43. 

As [NeV]A3425 limited the search to z > 0.109 we 
then scaled back the search to the three criteria excluding 
3. This search found one additional CLiF AGN candidate: 
J1241-f44. 

These search criteria failed to find the known CLiF 
AGN HI Zw 77, although it is in the SDSS database. We 
obtained the SDSS spectrum for HI Zw 77 separately 

While somewhat unsuccessful, this exercise serves to 
highlight both the apparent rarity of CLiF AGN in the AGN 
population, and the complexity of selecting such objects. 

In this paper we make use of SDSS spectra (Figure]^. 
The basic properties of the CLiF AGN sample studied in 
this paper are given in Table Note that the outputs of 
the SDSS pipeline are used only for the sample selection. 
Detailed measurements of emission line parameters such as 
the flux and velocity widths are measured using our own 
methods (@- The redshifts were determined using single 
Gaussian fits to the [OHIJA5007 emission line. This line was 
chosen because it is the most prominent emission line in the 
optical spectra of these and most other AGN. 


4 EMISSION SPECTRA 

In Figure we plot the spectra for the four CLiF AGN 
with SDSS spectra. There is a remarkable number of emis¬ 
sion lines. We performed a uniform line detection process by 
making fits to the lines and identifying them as described 
below. 

4.1 Line detection 

All line detections were determined by fitting single Gaus- 
sians to the emission features in the CLiF AGN spectra, us¬ 
ing the Starlink package dipscB accepting those with S/N > 
3. We dehne S/N as the total flux measured using the Gaus¬ 
sian prohle divided by the uncertainty of this flux from the 
fit as provided by DIPSO. For the three CLiF AGN without 
SDSS spectra (Q1131-I-16, ESO 138 G1 and Tololo 0109- 
373), we use the spectra and data presented in Table 
30-99 lines were found in the CLiF AGN (Table [^. 

® http://starlink.jach.hawaii.edu/starlink 
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III Zw 77 



Mrk 1388 



Figure 1. Labeled SDSS spectra for III Zw 77 and Mrk 1388. FHILs are labeled in red and marked with *. Note that we could not 
label all the emission lines on the spectra due to space constraints. The flux scale is measured in units of 10“^^ erg s“^ cm~^, the 

Observed Wavelength is measured in units of A. CONTINUED ON next PAGE. 
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J1241+44 



Observed Wavelength (Ang.) 

J1641+43 



-16 erg s-l A-l cm-2. 


Figure 1. Labeled SDSS spectra for J124H-44 and 1641+43. The flux scale is measured in units of 10 
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Table 2. Key results from the optical spectra of the CLiF AGN. Lines' refers to the total number of emission lines with S/N > 
3, ^FHIL is the number of detected FHILs, ^[FeVii] i® number of [FeVII] lines, Nfs is the total number of iron emission lines, 
[FeVII]/H^ is the flux ratio for the [FeVII]A6087 emission line, [FeX]/H^ is the flux ratio for the [FeX]A6375 emission line, 
[NeV]/H/3 is the H/3 flux ratio for the [NeV]A3425 emission line, Nj// is the number of detected HI lines and N? is the number of 
unidentified lines. The column ‘BL’ indicates whether broad emission lines were detected in the optical spectrum. ‘Y’ means yes, ‘?’ 
indicates tentative evidence and ‘N’ indicates a non-detection. All results were obtained using the Starlink package DIPSO. 


Name 

Lines 

^FHIL 

Ip > 54.4 eV 

N[FeVir] 

Nj.e 

[FeVII]/H/3 

[FeX]/H/3 

[NeV]/H/3 

Njrr 

N? 

BL 

Q1131-I-16 

99 

19 

7 

23 

0.50 

0.23 

1.01 

10 

30“ 

7 

ESO 138 G1 

71 

24 

7 

26 

0.24 

0.08 

1.22 

9 

2 

N 

Mrk 1388 

62 

13 

5 

11 

0.34 

0.19 

b 

7 

17 

7 

III Zw 77 

51 

13 

6 

11 

0.39 

0.21 

b 

8 

6 

Y“ 

J1241-I-44 

34 

12 

4 

9 

1.11 

1.88 

b 

3 

3 

N 

Tololo 0109-383 

31 

11 

5 

8 

0.28 

0.26 

0.57 

6 

d 

N 

J1641-I-43 

30 

8 

4 

7 

0.50 

0.21 

1.73 

4 

1 

N 


“ A few ‘unknown’ lines have been identified since [Rose et al.| ( |20TT] |. 

^ [NeV]A3425 is shifted out of the observed frame of the spectra for these objects. 

^ The BLR emission of III Zw 77 is significa ntly blueshifted (-675di25 k m s“^) with respect to the NLR emission. 
^ There were no ‘unknown’ lines indicated in [Fosbury &: Sansom| l [l983[ |. 
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Figure 2. Top. Multi-component model fits to the H/3 and 
[OIII]AA4959,5007 emission lines for 1241-1-44. The model fits are 
drawn with the dashed red line. The flux is measured in units of 
erg s“^ cm~^ and the wavelength is in units of A. Bottom. 
The residuals from the single Gaussian fits. 


Single-Gaussian profiles, did not always provide ade¬ 
quate fits to the wings or peaks of the strongest lines (EW 
> 10 Afl It is clear that these fits are inadequate for most 
of the CLiF AGN as the fitting residuals from the single 
Gaussian fits are large for the strong lines (Figures [3][^. 

To address this problem, double-Gaussian models were 
fitted to the strong spectral lines. To produce the models, 
both broad and narrow Gaussians were fitted to the H/3 
and [OIII] line profiles, where the narrow component widths 
for the H/3 and [OIII] lines are linked, yet the line cen¬ 
ters, intensities and broad H/3 width are free parameters. 
All but J1241-I-44 required a double-Gaussian model to fit 
these emission lines. Figures [2|5| show the residuals for the 
one and two-Gaussian fits to the [OIII]AA4959,5007 emis¬ 
sion lines. The variations shown in the residual plots of the 
fitting models used to obtain the emission line fluxes are 
dominated by noise in the continua of the spectra. Intrinsic 
velocity widths for the Gaussian components were obtained 
by correcting each of the components’ FWHM in quadrature 
using the SDSS instrumental width of ~3A. The intrinsic 
kinematics of the narrow (N) and broad (B) components 
are given in Table 

The parameters from the double-Gaussian fits were then 
used to fit all the strong (EW > 3 A) emission lines through¬ 
out the CLiF AGN spectra (FWHM, velocity shift and rel¬ 
ative amplitudes). 

The fits for the weakest lines (EW ^ 3 A) were overesti¬ 
mated with the double-Gaussian model because the broader 
component required by the double-Gaussian model tended 
to fit the local continuum. In these cases a single-Gaussian 
fit was sufficient. We indicate the number of Gaussians re¬ 
quired to fit each individual emission line in Tables |AHA4| 

We also take in to account constraints that are fixed by 
atomic physics. The [OIII]AA4959,5007 doublet was mod¬ 
eled using the constraint that they have the same EWHM, 
a 1:3 intensity ratio and the line center of A4959 fixed rel¬ 
ative to the center of A5007. The same approach was used 

^ These include Ha, H/3, [Oil]AA3726,3729 

[OIII]AAA4363,4959,5007, [NII]A6583, and [FeVII]A6087. 
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Figure 3. Top. Multi-component model fits to the H/3 and 
[OIII]AA4959,5007 emission lines for III Zw 77. The model fits 
are drawn with the dashed red line. The flux is measured in units 
of erg s~^ cm“^ and the wavelength is in units of A. Mid¬ 

dle. The residuals from the single Gaussian fits. Bottom. The 
residuals from the double Gaussian fits. 



4950 5000 5050 5100 5150 

Wavelength (Angstroms) 

Figure 4. Top. Multi-component model fits to the H/5 and 
[OIII]AA4959,5007 emission lines for Mrk 1388. The model fits 
are drawn with the dashed red line. The flux is measured in units 
of erg s“^ A~^ cm~^ and the wavelength is in units of A. Mid¬ 
dle. The residuals from the single Gaussian fits. Bottom. The 
residuals from the double Gaussian fits. 


for the [OI]AA6300,6364 (3:1 ratio) doublet blend (blended 
with [SIII]A6312 and [FeX]A6375) and the [Nil] AA6548,6584 
doublet (1:3 ratio), which is also often blended with the Ho; 
emission line. 

Throughout this paper, we use the summed fluxes of the 
two Gaussians. We do not use the data for weak emission 
lines (EW ^ 3 A) in any of the following investigation. The 
measured fluxes for the emission lines are presented for the 
SDSS CLiF AGN in Tables [ATPH in the Appendix. 


4.2 Line identification 

We regarded an identification to be secure for an emission 
line if its line center was within 1 cr of the wavelength pre¬ 


dicted for that particular emission line based on the redshift 
(as measured using the [OIIIJA5007 emission line). We used 


Osterbrock (19811 to identify the majority of the emission 


lines in the CLiF AGN spectra. For the rest of the peculiar 
emission lines we used Meinel, Aveni & Stockton (19751 and 


Kaler (19761. 


4.3 Emission line properties 


The emission lines are of a wide variety including 9-26 FHILs 
and 1-30 unidentified lines (Table [^. A full list of the mea¬ 
sured wavelengths, velocity widths (FWHM) and velocity 
shifts (Av) for the emission lines are given in Tables A1|A4 


in the Appendix. Below we discuss each group of emission 
lines. 
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Wavelength (Angstroms) 

Figure 5. Top. Multi-component model fits to the H/3 and 
[OIII]AA4959,5007 emission lines for 1641-1-43. The model fits are 
drawn with the dashed red line. The flux is measured in units of 
erg s“^ cm“^ and the wavelength is in units of A. Mid¬ 

dle. The residuals from the single Gaussian fits. Bottom. The 
residuals from the double Gaussian fits. 


4-3.1 Unidentified lines 

CLiF AGN host unknown emission lines in their spectra. 
Such emission lines are not reported in typical AGN. Four 
GLiF AGN have 1-3, but III Zw 77, Mrk 1388 and Q1131-I-16 
have 6, 17 and 30 unknown lines respectively (Tablej^. Only 
Q1131-I-16 and ESO 138 G1 share the same unknown emis¬ 
sion lines. Q1131-I-16 has 2 and 3 in common unknown emis¬ 
sion lines with J1241-I-44 and Mrk 1388 respectively. Mrk 
1388 has 1 and 3 in common unknown emission lines with 
J1241-I-44 and III Zw 77 respectively. Finally, III Zw 77 and 
J1641-I-43 share 1 common unknown emission line. 


Table 3. The double-Gaussian model parameters derived from 
the [01II]AA4959,5007 and H/3 emission lines of the GLiF AGN 
spectra. The ‘FWHMjv’ and ‘FWHMb’ columns present the rest- 
frame widths of the narrow and broad components of the model, 
which have been corrected for instrumental broadening. The ‘An’ 
column gives the rest-frame velocity separation of the broad rel¬ 
ative to the narrow kinematic components. The ‘Ratio’ column 
give the amplitude ratio for the intensities of the broad to narrow 
components. 


Name 


FWHMjv 
( km s“^) 


FWHMs 
(km s“^) 


Av Ratio, 
(km s“'^) 


HI Zw 77 
Mrk 1388 
J1241-I-44 
J1641-I-43 


216±11 

250±21 

145±9 

563±47 


512±32 

723±57 

1770±148 


-127±11 

-106±10 

-106±10 

-412±38 


0.29 

0.79 

0.30 


4-3.2 FHILs 

Each CLiE AGN has at least 10 FHILs. These FHILs 
are dominated by Fe ions, particularly [FeVII] (Ta¬ 
ble [2| Tables AllA4||. All CLiF AGN show emis¬ 
sion from [FeVII]AAAA3759,5159,5720,6087, [FeX]A6375 and 
[CaV]A5307. The [FeVII] A6087 line is particularly promi¬ 
nent, with flux ratios [FeVH]/H/3 of 0.24-1.11 (EW of 5-27 
A). However, in the spectra that cover [NeV]A3426, this line 
has flux ratios [NeV]/H(? of 0.57-1.73. 

In addition, J1241-I-44 has [ArV]A6435 and [FeXI]A6985 
that are not detected even in Q1131-I-16 and ESO 138 Gl, 
the two CLiE AGN with the most FHILs. 


4-3.3 Balmer emission lines 

Strong (EW > 3 A) higher level HI Balmer emission lines 
(He - Hie) are detected in 5/7 of the CLiF AGN (Table 
[^. Hydrogen emission from transitions with energy levels 
greater than He is not typically seen in AGN, implying some 
unusual physics in CLiF AGN, far from Case B. We inves¬ 
tigate the HI emission of CLiF AGN in more detail in §5.1| 


4-3-4 BLR emission lines 

In this work line width was not a selection criterion yet 
broad emission lines are not clearly detected in the spectra 
of the CLiF AGN suggesting a possible CLiF-type 2 AGN 
link. The exception is III Zw 77 in which the BLR Balmer 
emission lines are peculiarly blueshifted by -675±25 km s“^ 
with respect to the NLR Balmer emission lines. This is an 
expected signature of a black hole recoil event, where after 
the merger of two SMBHs, anisotropic emission of gravita¬ 
tional waves gives a ‘kick’ to the resulting SMBH after the 
merger (e.g. see Favata, Hughes fc Holz|2004 Civano et al. 
20121. While this phenomenon is interesting in its own right, 


III Zw 77 is the only recoil candidate to date that is a CLiF 
AGN. This subject is beyond the scope of this work. 


Rose et al. (20111 reported tentative evidence for broad 


emission at the base of the Hq:-|-[NH] blend for Q1131-bl6. 
In this work we assume that there is tentative evidence 
for a BLR if the Balmer recombination lines show broad 
wings that cannot be accounted for by using [OHIJA5007 
as an NLR template. Only Mrk 1388 shows such evidence. 
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The velocity width (FWHM) for this emission component 
is 1070±160 km s“^, however it is weak (represents 8% of 
the flux of the blend) compared with the narrow Ha+[NII] 
blend and is not present for higher order Balmer emission. 


5 RESULTS 


5.1 A unique CLiF region on the BPT diagram 

The Baldwin, Phillips & Terlevich (BPT) diagnostic dia¬ 
grams use emission line ratios to determine whether the 
emission region of an object is photoionized by an AGN, or 
by a starburst ||Baldwin, Phillips & Terlevich[1981| Veilleux 

- 4 


|fc Osterbrock|1987| Kewley et al.|2006 ). In Figure|6| we plot 

logic ([OIII]/H,5) vs logic ([Nil]/H q) for the CLiF AGN and 
SDSS emission line objects which include both AGN and 
starburst galaxies. The measurements for the SDSS objects 
were a subset (6000 objects) of the MPA-JHU SDSS DR7 
release of spectrum measurement^ This subset is composed 
of measurements from objects with spectra which were re¬ 
ported to have the highest S/N, as well as having [OIII]/H/3 
and [Nil]/Ha ratios which were consistent with the ranges 
shown in Figure]^ 

The SDSS AGN divide cleanly in the right-hand branch 
from star-forming galaxies in the left-hand branch. The solid 
line represents the extreme upper starburst line determined 
using the upper limit of theoretical pure stellar photoioniza¬ 
tion models ( Kewley et al.|20M I. The dashed line indicates 
the empirically derived boundary between HII region like 
galaxies and AGN ( Kauffmann et al.|2003 |. 

All the CLiF AGN fall in the cleft between AGN and the 
starburst galaxy branches. They are clearly removed from 
the main body of the SDSS AGN (grey dots) along the x- 
axis by logic [NII]/Ha ~ 0.2-1. Either the [Nil] emission 
is weaker, or Ha is enhanced in CLiF AGN compared to 
typical AGN. 

Balmer line recombination ratios are the most common 
tool to correct AGN spectra for intrinsic dust extinction as¬ 
suming conditions close to Case B. The Ha/H/3 ratio is the 
favored Balmer decrement for this task because these emis¬ 
sion lines are the strongest and have the widest wavelength 
separation. However, the large number of Balmer lines from 
high energy transitions shows that Case B does not apply 
to CLiF AGN ( Rose et al.|201l1 . 

In Figure we plot Hy/H/I versus Ha/H/J for the CLiF 
AGN (red circles) and the type 2 AGN of Reyes et al. (2008; 
black crosses). We limit the Reyes et al. (20081 sample to z < 


0.42 so that Ha lies in the SDSS spectrum, and use Balmer 
lines which have EW > 3 A, leaving 325 out of the original 
887 type 2 AGN. We use the SDSS pipeline measurements 
for the Balmer decrements of the type 2 AGN of [Reyes et| 
[aL]p008| ). 

A reddening line for [Calzetti, Kinney fc Storchi-| 


Bergmann (19941 dust is shown with Av= 0, 1 and 2 
marked. Most of the type 2 AGN of [Reyes et al.j ( [2008[ ) 
roughly follow the reddening trend indicated by the red¬ 
dening line. Within the la uncertainties 83% have values 
consistent with the reddening line. 

It is clear that while the Ha/H/3 ratios for the CLiF 


http://www. mpa-garching.mpg.de/SDSS/DR7 / 


01 


o 


QB 

O 



- 2-10 1 
Log([NII]/Ha) 


Figure 6. Diagnostic plot of Logio([OIII]/H/3) vs 
Logio([NII]/HQ!). AGN are defined to lie above the solid 
line, Hll-region like galaxies below the dashed line, and compos¬ 
ite galaxies between these boundaries. The red circles represent 
the CLiF AGN. The small points represent objects detected by 
the SDSS and are taken from the SDSS-DR 7 release. The region 
outlined by the red dotted line is a ‘CLiF AGN Region’. 



Figure 7. Plot showing H7/H/3 versus Hq/H/ 3. The CLiF AGN 
are indicated by red circles and the type 2 AGN of [Reyes et[ 
l [2008[ l are indicated by black crosses. The theoretical Case 
B values for the decrements are indicated by the dashed lines. 
The dark green squares indicate the Balmer decrement values for 
reddening of Ay = 1 and 2 (labelled). 


AGN all exceed the Case B values of 2.9 by up to a factor 
~2, the majority of the H 7 /H /3 agree with the Case B value 
of 0.47 within the uncertainties (Table 0. These Case B 
values apply to the low density (ue = 10 - 10 ^ cm“®) limit, 
with a temperature of lOOOOK ( jOsterbrock fc Ferland|2006 1 . 

Instead, 5/7 of the CLiF AGN are inconsistent with the 
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Table 4. The Hq/H/ 3 and H7/H/3 Balmer decrements for the 
CLiF AGN. 


Name 

Ha/H/3 

H7/H/3 

Q1131+16 

5.00±0.17 

0.44±0.02 

ESO 138 Gl 

4.25±0.10 

0.48±0.20 

Mrk 1388 

3.95±0.12 

0.49±0.03 

III Zw 77 

4.15±0.13 

0.34±0.04 

J1241-I-44 

6.36±0.21 

0.47±0.04 

Tololo 0109-383 

5.77±1.00 

0.48±0.05 

J1641-I-43 

5.38±0.19 

0.40±0.03 

Case B 

2.9 

0.47 


reddening line. These objects all have Ha/H/? ratios exceed¬ 
ing Case B, yet their H7/H/3 ratios are consistent with the 
Case B value. Only III Zw 77 and J1641-I-44 lie close to the 
reddening line and of these only III Zw 77 > 2a from case 
B H7/H/3. It is likely that reddening does not dominate the 
CLiF AGN Balmer decrements. 

A visual inspection of the spectra of the remaining 17% 


of the type 2 AGN of Reyes et al. 120081 that do not fol 


low the reddening line show that they are not CLiF AGN. 
Their positions on Figure are possibly due to uncertain¬ 
ties with the SDSS pipeline fits e.g. degeneracies with the 
H7-h[OIII]A4363 blend. 


5.2 Kinematics 

We investigated the kinematics of the CLiF AGN with SDSS 
spectroscopy using the velocity shifts (Av, and FWHM from 
the single Gaussian fits to the emission lines of the CLiF 
AGN because there are not enough lines for a broad and 


narrow component separation (14.ll. 


5.2.1 FWHM vs Ip 

In Figure]^ we plot the velocity width (km s“^) versus Ip 
(eV) for the emission lines detected in each of the CLiF 
AGN. There is significant scatter in the measured widths of 
a few hundred km s“^. However there is no trend for FWHM 
to increase with Ip. 

There are four lines (blue crosses) which have much 
higher FWHM than the median of all the emission lines of 
the CLiF AGN (horizontal dashed lines). These emission 
lines are: 

• [ArV]A7004 (Ip = 59.81 eV) in HI Zw 77 which may be 
intrinsically broader, but could always be blended with an 
unknown emission line. 

• [GaV]A5309 (Ip = 67.30 eV) in HI Zw 77 and J1241-h44 
possibly blended with [FeXIV]A5303. 

• H7 blended with [OH]A4415 in J1641-I-43. The low S/N 
of the blending companions leads to Gaussian fitting degen¬ 
eracies which are not resolved. 


The 3 GLiF AGN without SDSS spectra also show no 
trend for FWHM with Ip | Rose et al.|2011 Fosbury & San- 


|som||1983| and [Alloin et al.||1992p . However, [Alloin et al. 
|T992j~found a median velocity width for all the emission 
lines (FLILs and FHILs) around 1000 km s“^ for ESO 138 
Gl, 2-3 times higher than found for the SDSS CLiF AGN 
sample. 


5.2.2 Av vs Ip 

In Figure|^we plot the velocity shifts (Av; km s“^) relative 
to the [OIH]A5007 line, versus Ip (eV). Again there is signif¬ 
icant scatter of a few hundred km s“^. There is no evidence 
for a trend with Ip. 

There are 7 lines which have much larger Av (blue 
crosses) than the median of all of the emission lines of the 
CLiF AGN (horizontal dashed lines). These emission lines 
are: 

• [OH]AA3726,3729 doublet (Ip = 13.61 eV) in HI Zw 77, 
Mrk 1388 and J1241-I-44 whose blended line center (assumed 
to be 3727.4 A) is dictated by the electron density of the 
emission region. 

• [OI]A6364 (Ip = 0 eV) in Mrk 1388 which is in a blend 
with [FeX]A6375. 

• [FeXI]A6984 (Ip = 262.1 eV) in J1241-I-44 an uncom¬ 
mon [FeXI] emission line. 

• [FeVH]A5159 (Ip = 99.1 eV) in J1241+44 and 
[FeVI]A5176 (Ip = 75.0 eV) in J1641-h43 both of which are 
blended with each other. 


The outliers in blends are likely the result of fitting 
degeneracies. For the [FeXI]A6984 line, it is possible that 
this emission line has been misidentified. 


Q1131-I-16 also shows no trend for Av with Ip (Rose 


et al. 20111. Unfortunately Fosbury & Sansom (19831 and 


Alloin et al. (1992 1 do not comment on Av for the FHILs. 


5.3 The [OIII]A4363 Emitting Region 

The emission-line flux ratio of the FLILs [OHI] 4363/5007 
is a widely used diagnostic for the ionization mechanism 


and temperature of NLR in AGN (Osterbrock & Ferland 
|2006[ ) as it has little or no contamination from star formation 
compared to other NLR emission lines (e.g. [Nil], [SH] and 
[01]). However, [OHI]A43 63 has a 50 times higher cr itical 
density (nc = 10^ ® cm“®; Osterbrock & Ferland|2006 \ than 
[OHI]A5007 (nc = 10" ° 

FHILs. 


cm '^), comparable to those of the 


[Nagao et al. (20011 showed that, like [FeVH], the [OHI] 
4363/5007 ratio increases across Seyfert galaxy types from 
type 2 to type 1. This suggests that a significant portion 
of the [OHI]A4363 flux may be emitted from the [FeVH] 
emitting clouds in the CLiF region. 

We address this possibility directly in Figure where 
we plot the strength of [OHI]A4363 against [FeVH] A6087. 
Both values are normalized to [OHI]A5007 to eliminate AGN 
luminosity as a variable. Fi gure |10| includes both the GLiF 
AGN and the Erkens et al.| ([l9^ I FHIL sample. Erkens et 


al. (1997 \ observed a sample of 15 Seyfert galaxies and 2 
emission line radio galaxies based on the presence of FHILs 
in previous work. For both GLiF AGN, and the narrow-line 
[Erkens et al.| ( |1997[ ) samples, there is a clear correlation of 
[OHI]A4363 with [FeVH]A6087. The correlation coefficient 
for these objects is R = 0.92. Using Spearman’s rank corre¬ 
lation test, this correlation is statistically significant, with a 
p-value = 0.05%. 

This correlation clearly indicates that both the [EeVH] 
and [OHI]A4363 flrrxes are emitted from the same emission 
region in both GLiF AGN and in the narrow-line [Erkens et| 
[al.[ ( [1997[ ) AGN i.e. the GLiF region. Given that the FHILs 
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Figure 8. Velocity widths FWHM (km s”'^) for the single Gaussian fits to the emission lines versus Ip (eV) for the CLiF AGN (red 
circles). The blue crosses are for the emission lines which have a much higher FWHM than the median of all the emission lines of the 
CLiF AGN. The median velocity widths are indicated by the dashed lines. The dotted line indicates the boundary between the FLILs 
and FHILs (Ip = 54.4 ev; the Hell edge). 


in both the CLiF AGN and Erkens et al. (19971 narrow-line 


AGN samples share the same kinematics as the FLILs, this 
correlation suggests a strong link between the FHIL emitting 
clouds and the [OIII]A4363 emitting clouds. 

The broad-line [Erkens et al. (19971 AGN show no such 
correlation (R=-0.01, p=0.98). In the broad-line [Erkens etj 
al. (19971 objects FHILs do not share the same kinematics 


as the FLILs, but have higher FWHMs, and are blueshifted 
with respect to FLILs. The magnitude of the blueshift is cor¬ 
related with the FWHM. This is another difference between 
the CLiF AGN and typical type 1 AGN. 

The correlation of [OHIJA4363 and [FeVII]A6087 fluxes 
suggests that using the [OHI] 4363/5007 emission-line flrrx 
ratio is not a reliable way to estimate temperature, at least 
for the CLiF AGN. Whether this is a problem for typical 
AGN requires the determination of the level of contamina¬ 
tion of [OHIJA4363 from the GLiF region. 


5.4 Emission Region Diagnostics 

In order to investigate the physical conditions of both the 
FLIL and FHIL emission clouds, we u sed the photoion¬ 
ization code GLOUD\|^ ( [Ferland] [l998| ) to create plane- 
parallel, single-slab photoionization models for the CLiF 
AGN emission lines. We assumed that the clouds were radi¬ 
ation bounded and had a solar composition. 

We varied the ionization parameter (U) over the range 
-2.5 < logic [U] ^ 0, in increments of logic [U] = 0.5, and 
the hydrogen density (np) over the range 3.0 < logic [np 
cm“®] ^ 7.0 in steps of log [np cm“®] = 0.5. The models 
were run assuming ionizing continuum power-law indices a - 


The version of CLOUDY used in this paper is cl3.03. 
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Figure 9. Velocity shift (km s“l) with respect to the redshift of the CLiF AGN for the single Gaussian fits to the emission lines versus 
Ip (eV) for the CLiF AGN (red circles). The blue crosses are for the emission lines which have a much much larger Av than the median 
of all of the emission lines of the CLiF AGN. The median velocity shifts are indicated by the dashed lines. The dotted line indicates the 
boundary between the FLILs and FHILs (Ip = 54.4 ev; the Hell edge). 


1.5, -1.2 and -0.^^ We used the ‘table power law’ command 
of CLOUDY to describe the input ionizing continuum SED. 


This has a default spectral range from lO/rm to 50 keV (Fer- 
land||1998 b The results are shown in Figure [TI| We do not 
show grids which assumes a power-law of a=-1.5 or -0.8 be¬ 
cause these models failed to reproduce the observed emission 
line ratios for the CLiF AGN. 

We did not include dust grains when modeling the FHIL 
and FLIL clouds. This assumption is supported by the work 
of [Ferguson, Korista fc Ferland| 1 1997p,[Nagao et ah] ( |2003[ ) 


and Nagao, Maiolino & Marconi (20061 who demonstrated 


that for a large sample of AGNs the observed FHIL strengths 


The emission is described by F,, oc i/+“, where a is the power 
law spectral index. 


could not be reproduced in CLOUDY models that contain 
dust grains. 


5 . 4.1 FHILs 


11 


i) 


First we consider the FHILs. In Figure 
plot a grid of [FeVI] A5176/[FeVH]A6086 versus [Fe- 
VH]6086/5159 contours (iso-np and iso-U). We chose the 
[FeVI]A5176/[FeVH] A6086 ratio because it is sensitive to 
both U and the ionizing power-law, and we chose [Fe- 
VH]6086/5159 because this is the only FHIL ratio whose 
emission lines were detected in all CLiF AGN, that is sensi¬ 
tive to density yet not as sensitive to variations in U. Note 
that [FeVI] A5176 was not detected in HI Zw 77 and therefore 
it is not included on Figure [T^ a). 

In Figure [TI|[b) we plot a grid of [FeVH]3760/6086 ver- 
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Figure 11. (a) Diagnostic plot for the FLIL region showing the measured [FeVI]A5176/[FeVII]A6087 versus [FeVII]6087/5159 ratios for 
the CLiF AGN (red circles). Ill Zw 77 is indicated separately with a blue square. Iso density lines are solid and iso ionization parameter 
lines are dotted, (b). Diagnostic plot for the FLIL region showing the measured [FeVII]3760/6087 versus [FeVII]6087/5159 ratios, (c) 
Diagnostic plot for the FLIL region showing the measured [OI]A6300/[SII]A6717 versus [SIIJ6717/6731 ratios, (d) Diagnostic plot for the 
FLIL region showing the measured [NIIJA6583/H/3 versus [SIIJ6717/6731 ratios. 


sus [FeVII]6086/5159 contours. We chose [FeVII]3760/6086 
because it is sensitive to U. 

All the ratios studied in Figures |ll[ a)fc(b) consist of 
Fe lines. Therefore, they are completely free from the uncer¬ 
tainty in both the relative elemental abundance ratio, and 
the dust depletion factor. 

Comparing the line ratios for the CLiF AGN with the 
grid models in Figures [TT| ( a) fc(b) we find the ratios are con¬ 
sistent with uh in the range 4.5 < logio (n/r cm“®) < 6.5 
and U in the range -2.0 < log [U] < 0. 

We cannot fit a t/ value for the FHIL emitting clouds of 
III Zw 77 (blue square). This is true for all assumed power- 
laws for the ionizing continuum. 


5.4-2 Ha emitting region 


The observed Ha/H/i ratio is very large (3.95-6.36 with the 
Case B value of 2.9) for a given Hq/H/l ratio in CLiF ACN. 
A plausible explanation is that the high density of the CLiF 
region in these objects, coupled with an extended partially 


ionized zone (Gaskell & Ferland 1984 Rose et al. 20111, 


will lead to significant collisional excitation of Ha in the 
partially ionized zones of clouds in the CLiF region. The 
other Balmer series emission lines, will not be enhanced as 
they have lower cross-sections for collisional excitation, as 


well as higher excitation temperatures (Gaskell & Ferland 


1984 Osterbrock fc Ferland|2006 Rose et al.|2011 1. 

The conditions needed to produce the high observed 
Ha/H/1 flux ratios were investigated using GLOUDY. 
In Figure 12 we plot a grid of Ha/H/3 versus 
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Log,„([FeVIIj6087/[0IIIj5007) 



Figure 10. The strength of [OIIIJA4363 versus [FeVII]A6087, nor¬ 
malized to [OIIIJA5007. The CLiF AGN are indicated by red cir¬ 
cles, the narrow-li ne|Erkens et ah ||1997[ ) objects by green squares, 
and the broad-line |Erkens et al.f j 1997[ l objects by open green tri¬ 
angles. The solid line represents the best-fit slope to the CLiE 
AGN. 


Figure 12. Diagnostic plot for the CLiF region showing the mea¬ 
sured Ha/H/3 versus [FeX]A6375/[FeVII]A6087 ratios for the CLiF 
AGN (red circles). Ill Zw 77 is indicated separately with a blue 
square. Iso density lines are solid and iso ionization parameter 
lines are dotted. 


[FeX] A6375/[FeVII] A6087 contours (iso-nn and iso-U). We 
chose [FeX]A6375/[FeVII]A6087 because this ratio is sensi¬ 
tive to both U and the ionizing power-law, and because 
both [FeVII]A6087 and [FeX]A6375 have relatively high ncr« 
(10^ ® and 10®'^ cni“® respectively) and are therefore not so 
sensitive to variations in the density. 

Comparing the line ratios with the grid model we find 
that the Ha/H/3 is consistent with densities in the range 6.0 
< logic (nj7 cm“®) < 7.5. These densities are higher than 


found in typical NLR (log (nn cm < 4.0; Osterbrock & 
Ferland||2006 1 , and 0.5-1 dex higher than the FHIL clouds. 


A high Ha/H)? ratios could also be produced by lower 
metal abundances and/or a harder ionizing continuum 
(Gaskell & Ferland 1984| Rose et al.||201]^ . However, ex¬ 
tremely low metallicities (Z < 10“^“ Zq) are required in 
order to reproduce these ratios. At these metallicities it is 
questionable whether strong Fe emission lines would be ob¬ 
served. In addition, altering the slope of the ionizing con¬ 
tinuum results in the models being unable to reproduce the 
observed flux ratios of the FHILs. 


The [FeX]A6375/[FeVII] A6087 ratios imply U in the 
range -1.0 < log \U] < -0.5. This falls in to the range in 
U found for the FHILs in j |5XT| 

The ranges for U and xih for the FHIL and Ha emitting 
clouds of each CLiF AGN are given in Table 


5 . 4.3 FLILs 

Next we consider the FLILs, i.e. the normal NLR lines. In 
Figure [T^c) we plot a grid of [OI]A6300/[SH]A6317 versus 
[SH]6716/6731 contours. We chose [OI]A6300/[SH]A6317 be¬ 
cause this ratio is sensitive to both the ionizing power-law 
and U, and we chose [SHJ6716/6731 because it is sensitive 
to Ue | |Osterbrock fc Ferland|[2006[ ). 

In Figure I ll[d) we plot a grid of [Nil] A6583/H/3 versus 


Table 6. Logic U and logic ujr ranges for the FLIL regions of 
the GLiF AGN. Both lower and upper subscripts indicate the 
lower and upper limits of the ranges respectively. 


Name 

"Slower 

upper 

^H,lower 

,upper 

Q1131+16 

-2.0 

0 

3.5 

4.5 

III Zw 77 

-1.5 

-1.0 

3.0 

4.0 

Mrk 1388 

-2.5 

-1.5 

3.5 

4.5 

ESO 138 G1 

-2.5 

-1.0 

3.0 

3.5 

Tololo 0109-383 

-1.5 

-1.0 

3.0 

4.0 

J1241-I-44 

-2.5 

-1.0 

4.0 

4.5 

J1641-I-43 

-2.5 

0 

3.5 

4.5 


[SH]6716/6731 contours. We chose [NHJA6583/H/I because 
this ratio is sensitive to both the ionizing power-law and 
U. The [NH]A6583 emission is sensitive to the metallicity 


I Storchi-Bergmann|1991||Storchi-Bergmann et al. 19981. De¬ 
creasing the metallicity will increase the required U needed 
to reproduce the [NHJA6583/H/? ratio. However this is a 
prominent FLIL ratio detected in all the CLiF AGN spectra. 

The ratios for the CLiF AGN (red circles) are consistent 
for riH in the range 3.0 < logic (njy cm~®) < 4.5 and U -2.5 
< log [[/] < 0. The ranges for U and njr for the FLIL region 
of each GLiF AGN are given in Table These densities are 
consistent with those found for the NLRs of typical AGNs 


(e.g. Osterbrock & Ferland 20061. The lower end of the U 


range is consistent with those found for the NLR of AGN 
(-2.5 < log [U] < -2.0; |Bennert|[2005| |. The FHIL densities 
exceed those found for the FLILs by 2 dex. While there 
is some overlap for the U required by the FLIL and FHIL 
ratios, overall the U required by the FHILs is higher by 0.5 
dex assuming a solar metallicity. 
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Table 5. Logio U and logic ranges for the FHIL and Ha emitting clouds of the CLiF AGN. Both lower and upper subscripts 
indicate the lower and upper limits of the ranges respectively. 


Name 

Slower 

FHIL 

^ upper 

,lower 

,upper 

"Slower 

Ha 

U tipper 

,lower 

,upper 

Q1131-I-16 

-1.5 

0 

5.5 

6.5 

-1.0 

-0.5 

6.5 

7.0 

III Zw 77 

- 

- 

- 

- 

-1.0 

-0.5 

6.0 

6.5 

Mrk 1388 

-2.0 

-1.5 

4.5 

5.5 

-1.0 

-0.5 

6.0 

6.5 

ESO 138 Gl 

-1.5 

0 

4.5 

5.5 

-1.0 

-0.5 

6.5 

7.0 

Tololo 0109-383 

-1.5 

0 

5.5 

6.5 

-1.0 

-0.5 

6.5 

7.5 

J1241-I-44 

-2.0 

-1.0 

5.5 

6.5 

-1.0 

-0.5 

6.5 

7.0 

J1641-I-43 

-2.0 

0 

6.0 

6.5 

-1.0 

-0.5 

6.5 

7.0 


5.4-4 Comparisons of FHIL and FLIL region conditions 

Overall, the range of both U and riH increase from the FLIL 
emitting region (-2.5 < log \U] < 0, 3.0 < logic (n/r cm“®) 
< 4.5) to the FHIL (-2.0 < log \U] < 0, 4.5 < logic (nrr 
cm“®) < 6.5) and boosted Ha emitting clouds (-1.0 < log 
[17] < -0.5, 6.0 < logic (nrr cm“®) < 7.5) for all the CLiF 
AGN. The median increase in U from the FLIL to the FHIL 
and boosted Hq emitting clouds is ~ 0.5 dex for the lower 
end of the U range, and ~ 1.0 dex for the upper end. The 
median increase for n^r is ~2 dex at lower end of the uh 
range, and 3 dex for the upper end (see Tables [|&§. 

In Figure we map the results of |5.4.H5.4.3| on to 
logic U versus logic nrr space. The lines connect the FLIL 
(black circles), FHIL (red asterisks) and Ha (green squares) 
values for individual objects. We use the median U and nn 
values to mark the positions of the CLiF AGN emission 
regions on Figure |13| the error bars give the full ranges 
for these values. Ill Zw 77 and Tololo 0109-383 have the 
same median U and n^ values for the FLIL emitting re¬ 
gion and therefore overlap. Also ESO 138 Gl, J1241+44 and 
J1641-I-43 have the same median U and riH values for the 
FHIL emitting clouds. Finally, Q1131-I-16 and Tololo 0109- 
383 have the same median U and nn values for the boosted 
Ha emitting clouds. Both the median U and riH values in¬ 
crease from the FLIL emitting region to the FHIL emitting 
clouds and then to the Ha emitting clouds for all the CLiF 
AGN. However, we note that the range in U for the FHIL 
and boosted Ha emitting clouds overlap in all but two CLiF 
AGN: Mrk 1388 and J1241+44. 



Figure 13. A comparison of the CLOUDY model results for the 
FLIL, FHIL and Ha emitting clouds for the CLiF AGN. We plot 
logic U versus logic og. The results for the FLILs are indicated 
by black circles, the results for the FHILs are indicated by red 
asterisks and the results for the Ha emitting clouds are indicated 
by green squares. The error bars present the U and nu ranges. 
The results for individual objects are connected by solid lines. The 
results from the NLRs of the |Bennert| ( [2005t sample are indicated 
by open blue circles for the type 1 AGN, and open green triangles 
for the type 2 AGN. The expected range in density for the NLR 
is indicated by the red dotted lines ||Osterbrock &c Ferland|2006[l. 


For comparison, we use the nn and U measurements 
from Bennert ( 2005| ) for the NLR of normal AGN from spa¬ 
tially resolved VLT/FORSl and NTT/EMMI spectroscopy 
for 7 objects (5/7 type 1 AGN, open blue circles; 2/7 type 2 
AGN, open green triangles). This work shows that the range 
of riH (2.5 < logic n_f/ < 4.1) measured for these NLRs are 
consistent with the range found for the FLILs of 5/7 of the 
GLiF AGN, as well as the typical density of the NLR (2.5 < 
logic njr < 4.0; jOsterbrock fc Ferland||2006| ). However, the 
U values for the Bennert (20051 measurements are consis¬ 
tent with only the low values found for the GLiF AGN (-2.6 
< logic U < -2.0). The difference could be due to aperture 
effects, where we are sampling emission regions which are 
closer to the SMBH in the GLiF AGN using the SDSS spec¬ 
troscopy, when compared with the spatially resolved spec¬ 
troscopy of Bennert (2005). 


5.5 The CLiF region radial distance 


Given the computed U and uh for the FHIL and boosted 
Ha emitting clouds GLiF AGN from §5.4| we can calculate 
a range of distances (Rclif) to the GLiF regions given the 
ionizing luminosity (L/ojv): 


RcLiF = ( 


Lion 


AtvUuhc (ho) 


1/2 


( 1 ) 


where U is the ionization parameter, nn is the hydrogen 
density, {hi') is the mean ionizing photon energ^j^ and c is 
the speed of light. Note that we assume the CLiF emitting 


We assume a value of 56 eV for (hi'), which is based on an 
ionizing power law of -1.2, with the photon energy limits of 13.6 
eV and 5 keV. 
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l°gio (Kclif' P=) 


Figure 14. Logio Rsjjs ranges versus logic ^CLiF ranges for 
the CLiF AGN (red lines). The dashed black line indicates a one- 
to-one ratio of Rgi/s and RcLiF- 


region includes both the boosted Ha emitting clouds and the 
strong FHIL emitting clouds, as both these unusual emission 
phenomena are seen for all CLiF AGN. 

L/oiv is estima ted from using both the bolo- 

metric correction of [Heckman et (|2004| {Lbol= 3500 


L[o/jr i, with a variance of 0.38 dex), and the Elvis et al. 
( |1994| ) SED (Rbol ~ 3.1i0.2 Ljqjv). Given that L/ojv is 
proportional to and that we could not confidently 

correct L[oji/] for dust extinction, the computed radial dis¬ 
tances to the CLiF emitting region must be considered lower 
limits. 

Using this method we find a range of radial distances for 
the CLiF emitting region: 0.04 < KcLiF < 32.5 pc (Table 
[^. The uncertainties given were computed using the un¬ 
certainties for the [OHIJA5007 flux because they were the 
dominant source of uncertainty in the calculation. We could 
not calculate distance range for ESO 138 Gl, because the 
spectrum in Alloin et al. (19921 was not flux calibrated. In 
addition, the distance range calculations for HI Zw 77 relies 
solely on the CLOUDY grid presented in Figure [T^ 

The lower end of this distance range is consistent with 
the inner torus MAGNUM reverberation results of lYoshii eti 


al. (20141, who find AGN-torus distances of up to 0.1 pc. 


The higher end is consistent with disk radii measurements 
based on IR interferometry and high-resolution molecular 
line observations (e.g. Ghou et al. 2007 Tristram et al. 20071, 
as well as results from fitting IR SEDs and spectroscopy 
I Alonso-Herrero et al. [2011 1. 

We cannot make a direct comparison with the expected 
BLR radius for the GLiF AGN because we do not detect the 
BLR in the CLiF AGN spectra. 

We also determine dust sublimation distances (Rst/s) 
for the GLiF AGN using the relationship given by |Elitzur| 
1 20081. We assume dust sublimation temperatures in the 
range 1200 < Tsub < 1800 K, consistent with the expected 
range for torus dust grains ( jBarvairds ] |1987| [Rodrfgu^ 
Ardila fc Mazzalay|2006 1. The uncertainties were computed 


using both the uncertainty on the bolometric corrections and 
the uncertainties for the [OHI]A5007 flux. 

Figure [T4| plots the logic Rst/s ranges versus the logic 
RciiF ranges (pc) for the GLiF AGN (red lines) given in 
Table The x-axis error bars represent the Rcl/f distance 
ranges. These ranges intersect with the Rsc/s ranges at the 
median RciiF value. The y-axis error bars represent the 
Rst/s distance ranges. The lowest Rsfs values are for Rst/s 
at a sublimation temperature of 1800 K, and the highest 
values are at a temperature of 1200 K. The Rsc/s ranges in¬ 
tersect with the KcLiF ranges at distances which assume 
a sublimation temperature of 1500 K. The dashed black 
line indicates the one-to-one relationship between Rcf/f 
and Rsc/s- Note that all Rcf/f ranges intersect with the 
one-to-one line. 

However there is a caveat, L^qiii] has not been cor¬ 
rected for aperture effects. Due to the limited coverage of the 
SDSS fiber diameter (3 arcseconds), a fraction of the total 
[OHI] flux will be lost. Therefore, Lioiii] could be stronger 
than observed, which will subsequently increase both Ljon 
and Lbol, and therefore increase both Rcuf and Rsss- 
But the increase in Rcs/f and Rsc/s will be the same, and 
therefore the comparison of both parameters remain valid. 

We find a good agreement between the ranges for Rcs/f 
and Rsss as both ranges overlap for the CLiF AGN. In ad¬ 
dition the RciiF ranges encompass the one-to-one ratios of 
Rsss and Rclzf- Although we use L[o//i] to calculate both 
RcsiF and Rsc/s, the ratios between these two parameters 
do not have to agree. 


6 DISCUSSION 

6.1 The CLiF Region - the Inner Wall of the 
Torus? 


Given that both the dust grain sublimation and GLiF emit¬ 
ting region distance ranges are consistent ()5.5l, it is likely 
that the CLiF emitting region is closely associated with 
the inner wall of the torus, as suggested by [Murayama fc| 


Taniguchi (1998a I and Rose et al. (20111. 


This CLiF location is supported by the observation that 
FHIL emission is rich in iron lines particularly [FeVH], which 
can be enhanced relative to the dusty NLR (|Mor, Netzer &: 


Elitzur|2009 1 by the evaporation of dust grains in the inner 


torus wall, releasing the iron locked up in the grains (Pier 
jfc Voit|1995l. 


Rose et ar](|201l|) suggested that Q1131-I-16 is observed 


with our line of sight at a relatively large angle to the torus 
axis, so we have a maximal view of the far wall of the torus, 
while the AGN continuum is still obscured by the near side 
of the torus. 

We emphasize that such a geometry is consistent with 
the relatively narrow FWHM and small velocity shifts of the 
FHILs, since the range of circular velocities associated with 
the torus will be restricted largely to transverse motions 
in the sky plane, and any vertical gas motions would be 
directed close to perpendicular to the line of sight ([Rose etj 
al.|2011| . 


































































Coronal-Line Forest AGN 17 


Table 7. Radial distance (RcLsf) and sublimation distance (Rsub) ranges for the CLiF AGN. 


Name 

^CLiF 

lower 

RcLzF 

upper 

R'Stib 

ISOOK 

15007^: 

R'sub 

1200K 

Q1131-I-16 

0.52±0.01 

16.49±0.38 

0.85±0.05 

1.51±0.09 

4.34±0.25 

III Zw 77 

0.41±0.01 

2.33±0.06 

0.68±0.04 

1.09±0.06 

1.94±0.11 

Mrk 1388 

0.33±0.01 

32.52±0.74 

0.16±0.01 

0.53±0.03 

1.52±0.09 

Tololo 0109-383 

0.04±0.01 

2.03±0.18 

0.19±0.01 

0.53±0.02 

2.38±0.14 

J1241-I-44 

0.09±0.01 

1.98±0.13 

O.lOiO.Ol 

0.16±0.02 

0.29±0.02 

J1641-I-43 

0.51±0.02 

28.65±1.29 

1.48±0.09 

2.43±0.14 

4.24±0.24 


6.2 Collisional excitation of Ha 


In §5.1| we showed that the CLiF AGN occupy their own 
region on the BPT diagram due to a factor ~1.4-2.2 boosted 
Hq emission. 

This boosting may be due to collisional excitation in 
the high density CLiF region coupled with an extended par¬ 
tially ionized zone, so that Hydrogen-like ions can be ionized 
| |Gaskell fc Ferland|1984||Rose et al.|2011| . The other Balmer 
series emission lines, are not enhanced as they have lower 
cross-sections for collisional excitation, as well as higher exci¬ 
tation temperatures ([Gaskell fc Ferland|1984[ [Qsterbrock fc| 


Ferland|2006| Rose et al.|2011 1 . To explain the high Hq/H/3 

ratio measured in the spectra of the CLiF AGN, requires 
a major (>35%) contribution to the line flux of Ha from 
the CLiF region. Indeed, in Figure [l2| we investigated this 
possibility using CLOUDY modeling and found a required 
density range of 6.0 < logic (n/r cm“®) < 7.5, which is higher 
than expected for the NLR (< 10^'® cm“®) and the FHIL 
emitting clouds (10^ ® < uh < 10®'® cm“®). 

As there is a difference between the densities of the 
boosted Ha and FHIL emitting clouds, but no major dif¬ 


ference in U (15.4.41, we suggest that both sets of emission 
clouds are located in the CLiF region (possibly the inner 


torus wall, see 16.11, but there is a range of densities in the 


CLiF region, where the highest density clouds produce the 
boosted Ha flux. 

In Figure]^ we defined a box (red dotted line) to indi¬ 
cate the ‘CLiF AGN Region’ of the BPT diagram (limits: 
-1.0 < logic ([Nil]/Ha) < -0.4, and 0.65 < logic ( [OHI]/H,g) 
< 1.1). Within these limits, the BOSS catalog ( |Ahn et al.| 
20121 has 1284 CLiF AGN candidates. This population is 


just 0.5% of the 259008 objects in the AGN region of the 
BPT diagram (limits: -0.4 < logic ([Nil]/Ha) < 1, and -0.1 
< logic([OIII]/H/3) < 1.5). 

This region is based solely on the CLiF AGN studied 
in this paper. It is possible that higher levels of collisional 
excitation of Ha could push logic ([Nil]/Ha) values lower 
than -1.0. These CLiF AGN candidates will be investigated 
in a further paper. 


7 CONCLUSION 

When compared to ‘typical’ AGN, the optical spectra of 
CLiF AGN show unusual properties in the sense that they 
are rich in emission lines, and have a remarkable variety of 
species. 

Our analysis of seven CLiF AGN has shown that: 

1. The Ha Balmer emission line is boosted by a factor 


~1.4-2.2 when compared to higher energy Balmer transi¬ 
tions, putting them in a distinct location on the BPT dia¬ 
gram. 

2. There is no difference between the velocity widths of 
the FLILs and FHILs. 

3. There is no evidence that the FHILs are blueshifted 
with respect to the FLILs. 

4. The [OIII]A4363 emission line is emitted from the same 
region as the FHILs in CLiF AGN as its emission flux is 
correlated with that of [FeVH]A6087. 

5. The FHILs of CLiF AGN are emitted from dense clouds 
(lO'*'® < uh < 10®'® cm“®) when compared to the NLR. 

6 . The boosted Ha emission lines of CLiF AGN are emit¬ 
ted from denser clouds (10®'® < uh < 10^'® cm“®) when 
compared to the FHIL emitting clouds. 

7. We calculate small radial distances (0.04 < Kclif < 
32.5 pc) to the CLiF emitting region, which are comparable 
to the sublimation distance for dust grains (0.10 < Rsc/s < 
4.3 pc). 

A plausible model is that these objects are being viewed 
with a specific geometry which gives a maximal view of the 
inner wall on the far side of the torus, not at such a large 
angle from the disk plane that the BLR and continuum are 
exposed. In this scenario the inner torus wall is the FHIL 
and boosted Ha emitting region, which consists of gas with 
a range of densities that are needed to produce both strong 
FHILs and boosted Ha emission. 

We plan to follow up the CLiF AGN by (a) obtaining 
both X-ray and UV observations for the CLiF AGN, (b) 
modeling both the strength of the emission lines and kine¬ 
matics of the CLiF emitting region, (c) using the BPT cleft 
to And more (Rose et al. 2015), and (d) modeling the geom¬ 
etry (Glidden et al. 2015). 
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Table Al. Line identifications for the SDSS spectrum of Mrk 1388. ‘Ip presents the ionization potentials for the securely identified 
emission species. Flux ratios are relative to and are not corrected for reddening. Note that the fluxes presented here are from the 
double-Gaussian fitting model. Line IDs, Av and FWHM are determined by fitting a single Gaussian to the emission lines. The total flux 
of the H/5 emission line is (2.65di0.19) x erg s~^ cm~^ A“^. The column ‘Gaussians’ indicates the number of Gaussians required 

to fit the profile of the emission line. 


Line ID 

Ip (eV) 

A (A) 

^SDSS (A) 

F/Ph/3 

EW (A) 

Av (km s 

FWHM (km s-i) 

Gaussians 

[FeVII] 

99.1 

3758.9 

3839.5±0.1 

0.16±0.02 

7.4 

-42.8it7.9 

337.9it20.1 

2 

HI 

0 

3770.6 

3852.1±0.5 

O.OSitO.Ol 

1.5 

3.0it43.0 

154.2it60.5 

1 

HI 

0 

3797.9 

3878.9±0.2 

0.06±0.01 

3 

-76.4itl7.0 

380.4it43.8 

1 

HI 

0 

3835.4 

3918.4±0.9 

0.04±0.01 

2 

10.7it73.3 

228.8itlll.6 

1 

[Nelll] 

41.0 

3868.8 

3951.9±0.0 

1.30±0.28 

74.2 

-33.0it3.3 

366.4it8.1 

2 

[FeV] 

54.8 

3891.3 

3974.7±0.3 

0.28±0.08 

16 

-47.8it22.6 

445.7it54.9 

2 

[Nelll] 

41.0 

3967.5 

4053.5±0.0 

0.44±0.10 

25.8 

25.1it2.7 

379.0it6.1 

2 

[SH] 

10.4 

4068.6 

4155.8±0.1 

0.09±0.02 

3.9 

-82.5itl0.2 

382.4it6.8 

2 

[SH] 

10.4 

4076.3 

4163.3±0.1 

0.11±0.03 

4.8 

-75.6itl2.5 

382.4it6.8 

2 

H<5 

13.6 

4101.7 

4189.3±0.1 

0.24±0.05 

10.5 

-71.4it25.1 

382.0it6.1 

2 

7 


- 

4310.5±0.5 

O.OSitO.Ol 

1.6 

- 

194.2it43.3 

1 

7 


- 

4322.0±0.3 

O.OSitO.Ol 

1.6 

- 

206.3it32.2 

1 

7 


- 

4335.5±0.4 

0.05it0.01 

2.6 

- 

376.7it55.4 

1 

[OH] 

13.6 

4317.2 

4411.5±1.0 

0.04it0.01 

2.1 

70.2it69.1 

300.5itl03.0 

1 

H7 

0 

4340.5 

4434.0±0.1 

0.49it0.03 

24.7 

-20.6it8.6 

439.9it26.7 

2 

[OHI] 

35.1 

4363.2 

4456.9±0.1 

0.56it0.10 

28.3 

-38.5it6.0 

352.3itl5.7 

2 

Hel 

0 

4471.4 

4569.4±0.6 

0.06it0.02 

2.9 

92.9it39.0 

492.2it79.8 

1 

7 


- 

4759.4±0.5 

0.04it0.01 

2 

- 

260.0it50.6 

1 

Hell 

24.6 

4685.7 

4786.4±0.1 

0.34it0.08 

17 

-30.2it3.4 

301.9it8.1 

2 

H/3 

0 

4861.3 

4966.0±0.1 

- 

30.8 

-17.9it3.8 

350.8it9.2 

2 

[OHI] 

35.1 

4958.9 

5065.9±0.2 

2.86it0.19 

88.2 

-4.5itl0.3 

341.5it24.9 

2 

[OHI] 

35.1 

5006.8 

5114.9±0.1 

8.57it0.21 

264.5 

- 

357.0it8.7 

2 

[FeVII] 

99.1 

5159.0 

5269.9±0.3 

0.09it0.02 

3 

-32.0itl7.9 

283.2it31.3 

2 

[FeVI] 

75.0 

5176.4 

5288.2±0.4 

0.18it0.05 

6.1 

-0.7it22.0 

94.4itl7.9 

2 

[Nil] 

14.5 

5199.5 

5311.5±0.4 

0.04it0.01 

1.5 

17.3itll.2 

243.5it38.3 

1 

[FeVII] 

99.1 

5276.4 

5390.0±0.5 

0.07it0.02 

3.5 

-18.8it27.6 

269.4it47.3 

2 

[CaV] 

67.3 

5309.1 

5422.8±0.4 

0.15it0.04 

7.5 

-56.3it24.4 

539.6it60.1 

2 

Hell 

24.6 

- 

5529.9±0.5 

0.06it0.02 

3 

- 

346.5it56.5 

1 

[FeVI] 

75.0 

5425.7 

5540.9±0.8 

0.06it0.02 

3 

-105.5it44.0 

464.6it90.2 

2 

7 


- 

5557.9±0.5 

O.OSitO.Ol 

1.6 

- 

250.3it41.2 

2 

7 


- 

5578.2±0.3 

O.OSitO.Ol 

1.6 

- 

109.4it23.4 

2 

7 


- 

5739.7±0.5 

O.OSitO.Ol 

1.6 

- 

213.2it58.9 

2 

[FeVI] 

75.0 

5677.0 

5799.3±0.7 

0.04it0.01 

1.5 

-14.3it39.4 

278.0it82.0 

1 

[FeVII] 

99.1 

5720.7 

5844.1±0.1 

0.20it0.03 

6.2 

-7.4it4.2 

285.2it9.5 

2 

[Nil] 

14.5 

5754.6 

5878.8±0.3 

0.09it0.02 

3 

-2.7itl7.8 

423.0it47.5 

2 

Hel 

0 

5875.6 

6002.3±0.1 

0.12it0.02 

3.6 

-12.8it6.0 

321.0itl4.7 

2 

7 


- 

6207.1±1.5 

0.05it0.01 

1.8 

- 

532.6itl59.1 

1 

[FeVII] 

99.1 

6086.9 

6217.7±0.1 

0.34it0.05 

11 

-30.4it4.5 

284.2it8.7 

2 

[OI] 

0 

6300.3 

6436.6±0.1 

0.54it0.09 

14.8 

9.4it4.2 

367.9itll.5 

2 

[SHI] 

23.3 

6312.1 

6447.8±0.4 

O.lOitO.03 

2.8 

-30.6itl8.5 

295.8it36.8 

2 

[OI] 

0 

6363.8 

6501.6±0.3 

0.15it0.03 

4.1 

16.7itl4.2 

293.1it32.6 

2 

[FeX] 

233.6 

6374.6 

6511.1±0.2 

0.19it0.03 

5.2 

-53.9itl0.3 

263.2itl9.4 

2 

[Nil] 

14.5 

6548.1 

6691.8±0.2 

0.49it0.08 

15 

100.1it7.0 

301.9it4.9 

2 

Ha 

0 

6562.8 

6704.2±0.0 

3.95it0.12 

121.1 

-16.6it0.7 

288.4itl.8 

2 

[Nil] 

14.5 

6583.4 

6726.4±0.0 

1.47it0.16 

45.1 

33.6it2.1 

299.9it4.9 

2 

7 


- 

6822.3±0.7 

0.04it0.01 

1.8 

- 

451.6it61.0 

1 

[SH] 

10.4 

6716.4 

6862.2±0.1 

0.23it0.04 

7.8 

29.9it3.8 

327.9it9.5 

2 

[SH] 

10.4 

6730.8 

6876.6±0.1 

0.37it0.09 

11.2 

19.4it2.7 

311.0it6.3 

2 

7 


- 

7064.8±1.1 

0.12it0.03 

2 

- 

319.2it86.6 

1 

7 


- 

7148.2±3.0 

0.13it0.04 

2 

- 

649.3itl82.9 

1 

[ArV] 

59.8 

7005.9 

7157.5±0.5 

0.13it0.04 

2 

9.5it23.5 

251.9it54.7 

1 

7 


- 

7197.5±1.2 

0.15it0.03 

2.1 

- 

624.5itl53.6 

1 

Hel 

0 

7065.7 

7217.9±0.3 

0.07it0.02 

2 

-18.6itl3.0 

294.1it31.7 

1 

[Arlll] 

27.6 

7135.8 

7290.1±0.1 

0.12it0.02 

4 

5.2it4.2 

355.0itl0.8 

2 

7 


- 

7309.4±0.8 

0.04it0.01 

1.5 

- 

428.4it82.4 

1 

7 


- 

7326.5±0.8 

0.04it0.01 

1.5 

- 

348.5it85.7 

1 

[OH] 

13.6 

7319.9 

7478.8±0.5 

0.15it0.03 

5.1 

31.2it20.0 

521.4it48.0 

2 

[OH] 

13.6 

7330.2 

7490.2±0.5 

0.18it0.06 

6.1 

68.7it20.1 

340.3it31.8 

2 

7 


- 

7537.2±1.4 

0.07it0.02 

2.5 

- 

454.1itl26.6 

1 

7 


- 

7773.2±0.5 

O.OSitO.Ol 

1.4 

- 

196.4it29.2 

1 

7 


- 

7919.9±1.2 

0.07it0.01 

2.3 

- 

677.1itl77.7 

1 
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Table A2. Line identifications for the SDSS spectrum of III Zw 77. ‘Ip presents the ionization potentials for the securely identified 
emission species. Flux ratios are relative to and are not corrected for reddening. Note that the fluxes presented here are from the 
double-Gaussian fitting model. Line IDs, Av and FWHM are determined by fitting a single Gaussian to the emission lines. The total flux 
of the H/3 emission line is (2.52±0.09) x erg s“^ cm“^ A“^.The column ‘Gaussians’ indicates the number of Gaussians required 

to fit the profile of the emission line. 


Line ID 

Ip (eV) 

A (A) 

^SDSS (A) 

F/Fh/3 

EW (A) 

Av 

FWHM 

Gaussians 

[Oil] 

13.6 

3727.4 

3854.0±0.1 

0.63±0.07 

8.7 

-28.2ib7.8 

413.3ibl5.4 

2 

[FeVII] 

99.1 

3758.9 

3887.3±0.1 

0.49±0.05 

6.8 

17.4ibll.l 

450.2ib27.5 

2 

[FeV] 

54.8 

3839.3 

3968.0±0.8 

0.07±0.01 

2.1 

-164.1ib61.2 

409.0ibl09.1 

1 

[Nelll] 

41.0 

3868.8 

4000.6±0.0 

1.33±0.24 

19.2 

-11.9ib2.9 

346.2ib6.9 

2 

[FeV] 

54.8 

3891.3 

4023.9±0.2 

0.41±0.08 

5.9 

-5.6ibl6.8 

486.2ib33.1 

2 

[Nelll] 

41.0 

3967.5 

4103.1±0.2 

0.16±0.02 

2.8 

21.3ibl2.2 

183.7ib26.7 

1 

7 


- 

4109.1±0.5 

0.49±0.12 

8.6 

- 

1388.5ib77.0 

2 

[SII] 

10.4 

4068.6 

4208.2±0.3 

0 .12±0.02 

2.1 

64.2ib21.7 

464.2ib53.3 

1 

[SII] 

10.4 

4076.3 

4216.2±0.3 

O.OSibO.Ol 

0.8 

63.9ib21.6 

463.2ib53.5 

1 

H5 

0 

4101.7 

4242.1±0.1 

0.25±0.05 

4.4 

35.9ib8.5 

381.2ib21.0 

2 


0 

4101.7 

4233.0±1.4 

0.59±0.03 

10.4 

-643.5ibl00.1 

2546.1ibl71.9 

2 

[Oil] 

13.6 

4317.2 

4465.9±0.5 

0.15±0.04 

3.5 

99.4ib32.6 

345.2ibl0.5 

2 

H7 

0 

4340.5 

4489.0±0.1 

0.34±0.04 

7.8 

34.7ib6.7 

343.0ibl0.4 

2 

H7b 

0 

4340.5 

4479.3±0.8 

0.70±0.05 

16.1 

-648.3ib55.4 

2328.9ibl68.4 

2 

[OIII] 

35.1 

4363.2 

4512.2±0.1 

0.45±0.05 

10.4 

13.9ib4.9 

340.9ibl0.3 

2 

Hel 

0 

4471.4 

4625.7±0.7 

0.04±0.01 

1.9 

119.3ib47.4 

391.5ib79.7 

1 

7 


- 

4818.0±0.5 

0.06±0.01 

1.6 

- 

279.6ib70.3 

1 

Hell 

24.6 

4685.7 

4845.8±0.1 

0.44±0.07 

7.0 

17.0ib8.2 

469.4ib24.4 

2 

[ArlV] 

40.7 

4712.6 

4902.0±1.2 

O.OSibO.Ol 

1.4 

8.4ib75.9 

113.5ibl31.9 

1 

H/3 

0 

4861.3 

5027.9±0.1 

- 

15.4 

50.3ib3.8 

386.3ibl0.4 

2 

H/3s 

0 

4861.3 

5016.6±0.7 

2.03±0.09 

31.3 

-624.8ib41.3 

2628.4ib87.8 

2 

[FeVII] 

99.1 

4893.4 

5061.5±0.8 

O.OSibO.Ol 

1.3 

70.9ib51.8 

163.9ib69.6 

1 

[FelV] 

30.7 

4903.1 

5070.0±1.6 

0.05ib0.01 

2.7 

-16.5ibl00.4 

651.9ib249.1 

1 

[OIII] 

35.1 

4958.9 

5128.0±0.1 

1.94ib0.19 

29.9 

1.3ib6.8 

289.3ib4.8 

2 

[OIII] 

35.1 

5006.8 

5177.5±0.0 

5.82ib0.29 

89.6 

- 

285.6ib4.7 

2 

[FeVII] 

99.1 

5159.0 

5334.9±0.4 

0.08ib0.01 

3.2 

-0.2ib20.9 

226.2ib35.0 

2 

[FeVII] 

99.1 

5276.4 

5456.9±0.7 

0.06ib0.01 

3.2 

29.2ib40.0 

308.4ib88.6 

2 

[CaV] 

67.3 

5309.1 

5489.0±0.9 

0.22ib0.05 

8.9 

-60.7ib52.7 

975.0ibl59.7 

2 

[FeVII] 

99.1 

5720.7 

5916.7±0.1 

0.23ib0.03 

5.6 

46.3ib7.5 

340.9ibl7.4 

2 

Hel 

0 

5875.6 

6077.7±1.2 

0.24ib0.04 

5.8 

86.4ib60.3 

697.3ib99.1 

2 

7 


- 

6281.8±0.7 

0.05ib0.01 

2.9 

- 

348.5ib8.8 

1 

[FeVII] 

99.1 

6086.6 

6294.9±0.1 

0.39ib0.05 

9.5 

21.0ib3.9 

347.6ib8.8 

2 

[OI] 

0 

6300.3 

6515.1±0.1 

O.lSibO.Ol 

3.0 

-2.2ib6.0 

236.8ibl2.2 

1 

[Sill] 

23.3 

6312.1 

6527.0±0.7 

O.OSibO.Ol 

1.0 

-15.7ib34.6 

187.6ib47.0 

1 

[OI] 

0 

6363.8 

6581.1±0.3 

0.04ib0.01 

1.1 

13.6ibl4.9 

188.3ib22.0 

1 

[FeX] 

233.6 

6374.6 

6592.8±0.2 

0.21ib0.05 

4.8 

36.7ib7.1 

415.8ibl8.8 

2 

[Nil] 

14.5 

6548.1 

6806.8±0.4 

0.28ib0.07 

5.3 

-50.4ibl8.0 

373.2ibl2.4 

2 

Ho 

0 

6562.8 

6787.9±0.1 

4.15ib0.13 

81.5 

58.7ib3.4 

374.4ibl2.5 

2 

Hqb 

0 

6562.8 

6773.0±1.0 

7.09ibl.84 

139.2 

-658.5ib46.7 

2517.4ib76.5 

2 

[Nil] 

14.5 

6583.4 

6770.3±0.4 

0.89ib0.21 

15.9 

-50.5ibl8.1 

375.6ibl2.5 

2 

[SII] 

10.4 

6716.4 

6945.7±0.1 

0.15ib0.03 

3.0 

11.7ib3.4 

189.9ib6.5 

2 

[SII] 

10.4 

6730.8 

6960.5±0.1 

0.17ib0.03 

3.1 

8.8ib3.3 

207.4ib6.3 

2 

[ArV] 

59.8 

7005.9 

7243.0±1.3 

0 .11ib0.02 

2.8 

-73.8ib55.3 

1090.2ibl59.2 

1 

7 


- 

7282.6±1.2 

0.12ib0.03 

3.0 

- 

990.5ibl35.5 

1 

Hel 

0 

7065.7 

7307.8±0.4 

0 .11ib0.02 

2.8 

49.5ibl7.9 

477.5ib47.8 

1 

[Arlll] 

27.6 

7135.8 

7379.3±0.2 

0.08ib0.02 

2.5 

8.5ib7.7 

170.6ibl2.3 

1 

[Oil] 

13.6 

7319.9 

7569.1±0.5 

0.06ib0.01 

2.4 

-16.7ib22.0 

298.9ib43.1 

1 

[Oil] 

13.6 

7330.2 

7581.1±0.7 

0.04ib0.01 

1.8 

36.7ib30.1 

234.9ib83.7 

1 

7 


- 

7628.1±1.4 

O.OSibO.Ol 

1.7 

- 

342.3ib80.1 

1 

[FeXI] 

262.1 

7891.8 

8161.5±0.3 

0.13ib0.04 

7.4 

21.7ibl2.1 

320.2ib28.9 

2 

7 


- 

8738.4±1.0 

0.05ib0.01 

2.4 

- 

304.2ib76.9 

1 
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Table A3. Line identifications for the SDSS spectrum of J124H-44. ‘Ip presents the ionization potentials for the securely identified 
emission species. Flux ratios are relative to and are not corrected for reddening. Note that the fluxes presented here are from the 
double-Gaussian fitting model. Line IDs, Av and FWHM are determined by fitting a single Gaussian to the emission lines. The total flux 
of the H/3 emission line is (4.96di0.30) x erg s“^ cm“^ A“^.The column ‘Gaussians’ indicates the number of Gaussians required 

to fit the profile of the emission line. 


Line ID 

Ip (eV) 

A (A) 

>^SDSS (F) 

P/Fhi3 

EW (A) 

Av 

FWHM 

Gaussians 

[Oil] 

13.6 

3727.4 

3885.9±0.5 

0.77±0.18 

7.7 

104.2±44.0 

433.1±92.7 

2 

[FeVII] 

99.1 

3758.9 

3917.7±0.2 

1.09±0.15 

13.9 

17.2±3.7 

282.4±40.2 

2 

[Nelll] 

41.0 

3868.8 

4031.7±0.1 

1.57±0.15 

13.7 

-27.8±8.8 

256.3±22.2 

2 

H 7 

0 

4340.5 

4523.2±0.3 

0.47±0.04 

3.4 

-31.8±18.6 

205.5±34.8 

2 

[OIII] 

35.1 

4363.2 

4547.1±0.2 

1.3±0.13 

8.3 

-17.2±10.8 

287.9±29.6 

2 

Hell 

24.6 

4685.7 

4883.4±0.1 

0.87±0.09 

4.5 

-3.1±9.3 

227.4±18.4 

2 

[Nelll] 

41.0 

4815.9 

5017.9±0.3 

0.13±0.04 

2.8 

-73.3±21.5 

89.7±38.8 

2 

H/3 

0 

4861.3 

5066.2±0.1 

- 

5.1 

-13.2±5.3 

178.5±12.8 

2 

[OIII] 

35.1 

4958.9 

5168.2±0.1 

1.39±0.11 

7.5 

2.4±4.8 

206.1±11.7 

2 

[OIII] 

35.1 

5006.8 

5218.0±0.0 

4.34±0.28 

22.6 

- 

215.8±4.7 

2 

[FeVII] 

99.1 

5159.0 

5375.0±0.4 

0.21±0.07 

2.1 

-97.7±22.4 

141.3±67.3 

2 

[FeVI] 

75.0 

5176.4 

5394.8±0.3 

0.79±0.18 

7.9 

-3.8±17.8 

147.0±38.9 

2 

[CaV] 

67.3 

5309.1 

5531.4±1.1 

0.27±0.07 

2.1 

-94.0±59.8 

564.4±136.9 

2 

7 


- 

5766.6±0.3 

0.29±0.06 

2.2 

- 

168.7±36.4 

2 

7 


- 

5856.0±0.4 

0.27±0.07 

2.1 

- 

157.2±47.3 

2 

[FeVII] 

99.1 

5720.7 

5963.0±0.2 

0.64±0.09 

4.1 

44.6±10.1 

184.4±27.7 

2 

Hel 

0 

5875.6 

6124.2±0.5 

0.17±0.04 

1.9 

30.6±12.1 

159.0±16.5 

2 

[FeVII] 

99.1 

6086.9 

6343.6±0.1 

1.11±0.08 

5.2 

-10.0±4.5 

209.2±11.0 

2 

7 


- 

6443.5±0.5 

0.11±0.04 

1.6 

- 

139.1±52.0 

2 

[OI] 

0 

6300.3 

6566.0±0.2 

0.61±0.06 

3.2 

-5.8±7.3 

186.9±17.5 

2 

[OI] 

0 

6363.8 

6632.4±0.3 

0.21±0.04 

1.3 

0.9±13.5 

147.9±2.4 

2 

[FeX] 

233.6 

6374.6 

6643.7±0.1 

1 .88±0.12 

8.7 

5.9±2.5 

219.9±6.2 

2 

[ArV] 

59.8 

6435.1 

6707.0±0.5 

0.13±0.04 

1.5 

16.0±25.4 

176.5±50.3 

2 

[Nil] 

14.5 

6548.1 

6823.7±0.7 

0.17±0.06 

1.3 

-29.9±30.0 

170.5±51.2 

2 

Ha 

0 

6562.8 

6839.7±0.0 

6.36±0.21 

30.5 

-3.1±1.1 

213.7±2.8 

2 

[Nil] 

14.5 

6583.4 

6861.2±0.2 

0.71±0.08 

3.9 

-2.9±7.8 

180.7±18.5 

2 

[SH] 

10.4 

6716.4 

6999.0±0.5 

0.11±0.03 

1 

-33.9±20.1 

78.1±97.2 

2 

[SH] 

10.4 

6730.8 

7014.8±0.4 

0.22±0.04 

2 

-3.3±17.8 

181.8±35.5 

2 

[FeXI] 

262.1 

6985.2 

7278.5±0.3 

0.24±0.04 

2.4 

-61.5±13.0 

142.2±25.1 

2 

[ArV] 

59.8 

7005.9 

7299.5±0.8 

0.24±0.06 

2.4 

-82.8±32.6 

260.6±65.5 

2 

[ArIH] 

27.6 

7135.8 

7436.6±0.4 

0.21±0.07 

2 

-13.3±16.6 

144.6±74.0 

2 

[OH] 

13.6 

7319.9 

7629.2±0.5 

0.32±0.08 

2.2 

14.1±20.3 

195.1±62.2 

2 

[OH] 

13.6 

7330.2 

7639.7±0.6 

0.21±0.06 

1.9 

6.8±25.5 

175.5±50.5 

2 

[FeXI] 

262.1 

7891.8 

8225.0±0.1 

1.45±0.13 

7.6 

4.3±5.6 

208.3±15.6 

2 
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Table A4. Line identifications for the SDSS spectrum of 1641+43. ‘Ip presents the ionization potentials for the securely identified 
emission species. Flux ratios are relative to and are not corrected for reddening. Note that the fluxes presented here are from the 
double-Gaussian fitting model. Line IDs, Av and FWHM are determined by fitting a single Gaussian to the emission lines. The total flux 
of the H/5 emission line is (2.26+0.10) x 10~^^ erg s~^ cm~^ A“^. The column ‘Gaussians’ indicates the number of Gaussians required 
to fit the profile of the emission line. 


Line ID 

Ip (eV) 

A (A) 

^SDSS (E) 

F/Fh/3 

EW (A) 

Av 

FWHM 

Gaussians 

[NeV] 

97.1 

3345.6 

4084.2±0.4 

0.46±0.18 

16.0 

-174.3±133.5 

744.5±99.6 

2 

[NeV] 

97.1 

3425.9 

4182.1±0.2 

1.73±0.13 

47.9 

-164.2±62.5 

858.8±37.1 

2 

[Oil] 

13.6 

3727.4 

4552.4±0.2 

0.48±0.02 

20.6 

-8.6±0.2 

490.3±41.3 

2 

[FeVII] 

99.1 

3758.9 

4588.8±0.8 

0.36±0.09 

15.2 

-149.5±222.3 

809.3±129.1 

2 

[Nelll] 

41.0 

3868.8 

4725.7±0.6 

0.95±0.04 

41.9 

22.1±1.3 

827±61.4 

2 

[FeV] 

54.8 

3891.3 

4752.2±0.8 

0.20±0.08 

8.3 

-40.0±3.8 

741±68.9 

2 

[Nelll] 

41.0 

3967.5 

4846.0±0.4 

0.30±0.04 

12.9 

6.9±0.3 

369.9±58.2 

2 

7 


- 

4854.0±0.8 

0.07±0.01 

3.1 

-102.3±111.9 

369.2±58.1 

1 

H5 

0 

4101.7 

5009.9±1.0 

0.18±0.04 

6.6 

7.2±0.7 

665.7±224.4 

2 

H7 

0 

4340.5 

5303.6±0.8 

0.40±0.03 

15.8 

122.0±168.0 

1222.0±124.4 

2 

[OIII] 

35.1 

4363.2 

5329.0±0.7 

0.39±0.04 

15.6 

-14.7±2.2 

819.2±77.7 

2 

[Oil] 

13.6 

4414.9 

5396.9±2.7 

0.15±0.04 

5.9 

251.5±120.0 

710.0±59.2 

2 

Hell 

24.6 

4685.7 

5721.2±0.7 

0.30±0.03 

9.6 

-102.3±111.9 

680.9±73.6 

2 

H/3 

0 

4861.3 

5937.8±0.3 

- 

35.0 

8.3±0.3 

721.3±39.8 

2 

[OIII] 

35.1 

4958.9 

6057.3±0.4 

1.77±0.08 

67.3 

26.9±4.3 

626.3±45.3 

2 

[OIII] 

35.1 

5006.8 

6115.4±0.1 

5.31±0.25 

201.9 

- 

692.3±16.5 

2 

[FeVII] 

99.1 

5159.0 

6302.5±0.2 

0.08±0.01 

3.5 

64.1±10.9 

546.0±24.8 

2 

[FeVI] 

75.0 

5176.4 

6315.9±1.1 

0 .11±0.02 

4.0 

-312.7±200.0 

543.4±104.7 

2 

[FeVII] 

99.1 

5720.7 

6985.9±0.9 

0.25±0.04 

10.2 

-59.3±63.6 

788.0±181.7 

2 

[FeVII] 

99.1 

6086.9 

7430.2±1.0 

0.50±0.12 

19.4 

-174.5±78.0 

904.1±92.9 

2 

[OI] 

0 

6300.3 

7698.9±2.4 

0.21±0.05 

6.6 

144.5±50.0 

648.2±237.2 

2 

[OI] 

0 

6363.8 

7778.2±3.4 

0.07±0.02 

1.7 

210.5±100.0 

954.5±316.7 

1 

[FeX] 

233.6 

6374.6 

7785.7±0.7 

0.21±0.05 

7.2 

-39.9±8.2 

814±85.2 

2 

[Nil] 

14.5 

6548.1 

7997.1±0.5 

0.57±0.1 

17.6 

-29.0±10.0 

364.4±69.9 

2 

Hq 

0 

6562.8 

8017.7±1.6 

5.38±0.19 

166.0 

60.8±15.8 

476.8±38.2 

2 

[Nil] 

14.5 

6583.4 

8040.2±0.5 

1.71±0.11 

52.8 

-29.4±10.3 

362.0±69.4 

2 

[SII] 

10.4 

6716.4 

8203.6±1.0 

0.19±0.05 

8.0 

5.8±0.7 

479.7±170.7 

2 

[SII] 

10.4 

6730.8 

8221.3±1.3 

0.30±0.04 

5.0 

10.6±3.2 

350.3±84.5 

2 


Table Bl. Emission lines without a secure ID for the SDSS spectra. Information for Q1131+16 is presented to highlight any overlap. 


^rest (A) 

Mrk 1388 

HI Zw 77 

J1241+44 

J1641+43 

Q1131+16 

3973.5 


Y 


Y 


4219.4 

Y 





4230.6 

Y 





4243.8 

Y 





4659.1 


Y 




5413.0 

Y 




Y 

5440.4 

Y 





5460.3 

Y 





5533.1 



Y 


Y 

5618.6 

Y 


Y 


Y 

6075.1 

Y 

Y 




6182.6 



Y 



6678.1 

Y 





6915.4 

Y 





6997.1 

Y 





7043.9 

Y 

Y 




7154.9 

Y 




Y 

7171.6 

Y 





7377.1 

Y 

Y 




7608.8 

Y 





7752.4 

Y 





8450.0 


Y 






